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ABSTRACT 


This thesis investigates the problem of aerial search for a mobile 
surface target. The initial position of the target is known, and a 
certain amount of time elapses before the search begins. The target 
has a variety of choices as to speed and direction, and the search air- 
craft has a choice of search patterns to be flown. The method proposed 
to solve the problem of maximizing or minimizing the probability of 
detection is the game theoretic analysis of various strategies avail- 
able to the searcher and the evader. Four strategies for each are 
compared and evaluated, with optimal strategies, within the set of 


strategies investigated, being derived. 
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La Sent roduction. 

A frequently encountered problem of attack aircraft is the detection 
and destruction of mobile enemy installations at some time after initial 
detection. For example, a flight may be fired upon by a mobile surface- 
to-air missile site. Upon return to the air base or carrier the location 
of the attack is reported and a strike against it is readied and sent out. 
In the meantime, the missile site has often had time to depart the posi- 
tion in which it was originally situated. There are many other ways in 
which a mobile target may be detected, such as aerial reconaissance, air- 
borne electronic detection, or visual sighting. Some of these may give 
a probabilistic rather than a deterministic position. 

The mobile target, if it knows it has been discovered, has the op- 
tions of remaining in its present poSition or moving. If it moves, it 
has some variety in its choices of direction and speed. Often these 
choices of movement are to some degree limited. The search aircraft are 
then confronted with the problem of finding and attacking the moving 
strike objective. This can be hampered by such variables as time to sun- 
set, terrain, weather, and other uncontrollable elements. The target may 
be found by methods other than visual if it is emitting electronic sig- 
nals, i.e., if it has not moved or if it has moved and is set up and 
operating in another position. The most difficult problem arises when 
the method of detection must be strictly visual, as when the objective 
is moving or if it has set up in another position and is not emitting. 
The degree of difficulty also depends on how much freedom of motion the 
target has. If the direction of travel is limited rather strictly, as 
when there are only a small number of roads along which it may travel 


and off-road travel is not possible, the searcher's task is somewhat 


simplified. If, however, it is free to travel in any or almost any direc- 
tion, the problem becomes much more difficult. This thesis examines some 
of the problems confronting a searcher faced with the problem of discover- 
ing a target with unrestricted direction of movement, and with varying 
maximum speeds. 

This problem may be attacked in two fashions. One method is to con- 
sider the movements of the target as being determined probabilistically, 
with such variables as direction, velocity, and time of movement all being 
described by some density function, so that an estimate of its most likely 
position, or most likely set of positions, can be found using the usual 
rules of probability. In fact, all these variables have some probabilis- 
tic limitations, but within these limitations there is another influence: 
the judgment of the commander of the site. It would be somewhat naive to 
think that the commander would be guided solely by the constraints imposed 
by the state of nature. Rather, it 1s more reasonable to assume that he 
would try to vary his behavior within the limits of his physical con- 
Straints so as to attain maximum 1iklihood of evasion. If this did not 
occur immediately, then certainly experience would lead him to adopt more 
sophisticated evasive.tactics. The intelligent searcher would then alter 
his search pattern to counter the evasive tactics to as great a degree as 
possible. This development brings the problem into the area of game 
theory, and it is with this aspect that this paper concerns itself. 

To simplify the problem a number of probabilistic variables are 
treated deterministically. The effect of weather, time of day, terrain 
and variance in terrain over the search area, and other external influ- 
ences are ignored. These and other assumptions are discussed in the next 


section. A small number of strategies for both searcher and evader are 


discussed. While an attempt has been made to consider the most advan- 
tageous strategies, there are probably many of merit which have not been 
included. It is felt that an approach has been made to the solution of 
the problem, however, and some possibly beneficial tactics derived. 

The general approach is to select several tactics, both for the 
searcher and the target. These tactics are then evaluated in pairwise 
combinations, using game theory criteria. This allows conclusions to be 
drawn concerning the relative merits of the strategies. 

The next section presents the baSic assumptions of the model and 
attempts to justify them. The ensuing sections develop the general 
theoretical framework and the specific mathematical formulations of the 
proposed strategies. After comparing these strategies, some conclusions 


as to the value of each are inferred. 


II. Theoretical Basis for the Model. 
2.1 Assumptions. 

It would be an almost impossible task to obtain a universal solution 
to this search problem, in which all relevant factors were included. On 
the other hand, elimination of a critical variable would leave the con- 
clusions worthless. The assumptions made in the development of the model 
are set forth here in an attempt to explain its limitations and clarify 
its objectives. 

The space in which the game is played is a two-dimensional Euclidean 
geographical space. The terrain problem is ignored in that the area is 
considered to be perfectly flat, and in no way affects either target mo- 
bility or the searcher's visibility. The target is free to move in any 
direction and at any velocity up to its acre i for any duration of time. 
The searcher has a fixed range sighting probability; within this range the 
probability of detection is 1.0, outside this range the probability of 
detection is 0. The searcher is capable of flying any pattern exactly. 
The initial position of the target is known exactly. The target knows 
he has been detected, and takes action accordingly. No variables such 
as weather, time of day, or other external phenomena are considered. The 
time at which the aircraft will arrive to begin his search is known to 
the target as soon as the initial detection occurs. 

There are other minor assumptions in the model. These cited are 
considered to be the most relevant. To be able to claim that the model 
has any relevance to the real world some of these assumptions must be 
justified. 

First, consider the playing area. Unless the terrain is extremely 


rough, the aircraft acts essentially as if it were flying above a two- 


10 


dimensional area. The behavior of the target is more affected by this as- 
sumption. Since the target-chosen factors are speed and direction, the 
effect of terrain on these is relevant. In a real-world situation, the 
target might indeed be somewhat terrain limited. However, most missile 
and radar sites are near cities, which in turn are usually situated in 
flat terrain. Also, there is usually a heavy density of roads near any 
city, which gives the target a great degree of flexibility in its move- 
ments, both as to direction and speed. If the area around the initial 
position is not travelable in almost any direction, the problem is simpli- 
fied to some extent, but is susceptible to the same type analysis as is 
proposed here. 

The assumption that the target can move at any speed up to its max- 
imum for the duration of the search is reasonable if the search aircraft 
arrive within a reasonable length of time, on the order of four hours 
after initial detection, and stay on station a moderate length of time. 
Since most attack circraft can fly at low altitudes for only about two 
hours, the total time the target must be able to keep moving is about 
six or seven hours. This is well within the capability of most mobile 
sites. 

The method of first detection of the target is often visual or photo 
graphic, in which cases the poSition is known exactly. If the position 
1s not known exactly, the point of highest probability is usually taken 
as the true position and treated deterministically in any actual search. 
The effect of large errors in the position of the original site is not 
examined in this discussion. 

The implications of a fixed range of detection are well known and 


no further implications have arisen in this analysis. {3] However, the 
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idea that an aircraft can fly any desired pattern exactly will strike any 
experienced aviator as a flagrant misconception. In fact, this is prob- 
ably the most critical assumption of the model. Recent years have seen 

the development of Doppler navigational equipment which enables an air- 
craft to fly a precise track, and airborne radar in a control aircraft 

can guide other aircraft over a desired pattern with good accuracy. There 
is no reason to believe that such techniques could not be applied to search 
alrecraft if the situation warranted. 

Another critical assumption is that the target knows, as soon as de- 
tected, at what time the search aircraft will arrive. Some examination 
of the actual procedures reveals this assumption to be less restrictive 
than it initially appears. If the target moves, the probability of de- 
tection diminishes with time in any of the strategies considered. The 
search aircraft will therefore want to begin the search as rapidly as 
possible. The target often knows the position of the base or carrier 
from which the attack aircraft will be launched, and can therefore pre- 
dict fairly accurately their arrival time. An estimate might also be ob- 
tained from observing past arrival times. 

While these assumptions restrict the application of the model to 
a limited situation, they allow the formulation of a general mathematical 
model. Refinements of the general approach may be made to deal with 
specific deviations. 

2.2 Outline of Game Theory Approach. 

The searcher, denoted player A, and target, denoted player B, have 
available to them a large and possibly infinite set of strategies. Each 
attempts to maximize his payoff by choosing judiciously among these 


Strategies. B may choose any direction and any speed up to his maximun, 
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and may stop his movement at any time. A may choose any Search pattern 
over the playing area which does not exceed his range. The payoff func- 
tion for A is the probability of target detection. The payoff to B is 
either the negative of the probability of detection, -P(D), or the prob- 
ability that A does not detect B, 1-P(D). If the former is defined as 
B's payoff the game is zero-sum; if the latter the game is constant-sum 
and the sum is one. 

The set of strategies for player A is denoted [’ and each element of 
this strategy set is denoted &. The set of strategies for player B is 
denoted 7 and the elements of the set are denoted 6 . The pairs of 
strategies (*%, 4) are elements of the Cartesian product of J” and a 
denoted fa x é . The payoff to A for a given pair of strategies, “ce 
and hs, , is MCX, A ) = P(A detects B/ X* %, 9 Pp = f8,)- For an opti- 


mal strategy to exist, there must exist strategies 2 and Ge such that 


min max M(&<,f) = max min M(%,f) = MC oR”, fp) (4) 
po - 
Whether in fact such strategies exist has not been determined. To 
determine the existance of such strategies, either all possible strate- 
gies must be examined, or it must be proven topologically that such strat- 
egies exist. The first method is difficult due to the large number of 
possible strategies. There exist some existence theorems which pertain 
to a pursuit-and-evasion game, but none were found which were applicable 
to the game under discussion. An existence theorem developed by Ryll- 
Nardzewski [1] which was the most general as to existence criteria, indi- 
cated that a solution involving pure strategies did not exist. 
It 1s also possible that an optimal solution which is a convex com- 


bination of various strategies exists. To find this, it is necessary to 


ils: 


examine the payoff matrix of the various strategies, and to apply the 
usual methods of game theory to find a solution. This requires develop- 
ing the probability of detection for various strategies. Therefore it is 
worthwhile to develop payoff functions for individually considered strate- 
gies of A and B. 

The general procedure then is to choose some strategy &%, and some 
strategy Po : 3, specifies a rule for choosing direction, speed and 


stopping point for the target. c¢, specifies a pattern to be flown over 


6 
the playing area. The payoff is P(A detects b/h=aA, 18 = 3, ). These 
strategies may be probabilistic and specify a probability function for 
choosing speed or initial aircraft search point, or may be deterministic, 
such as always start the search at a radius of five miles, or always re- 
main at the initial position. If the set of attractive strategies is 


compared pairwise in this manner, the matrix of strategy payoffs is gen- 


erated. 
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III. Development of Strategies. 
3.1 Strategies Considered. 

As has been noted, there are a great many possible strategies for 
both searcher and target. After considering a large number of strategies, 
the following were selected for investigation on the basis of their intui- 
tive attractiveness or because they are frequently used in practice. A 
glossary of terms is given in table one. The strategies as proposed by 


the author are: 


x 
iL 


Player A searches over Me in the most efficient manner which also 
searches over all directions of possible target travel, (0, 2). This is 
taken to be a spiral search, and the choice of the target velocities to 
be searched is determined in a game theoretic manner. 

Table I. Glossary of Terms 
ee target speed 

“a search aircraft speed 


Vv maximum target speed 


t 
max 


© direction of target movement or polar angular coordinate 
of aircraft position 

R aircraft range 

Tae minimum time after sighting at which target may move 

ah time from initial detection of aircraft arrival 


arr 
t <a - T , or time the target has been in motion at aircraft 
ae min 
arrival time 


r radius from initial target position 


To radius at which the aircraft Starts its search 


1 


th 


S; circumference of the 1 circuit searched by the aircraft 


W fixed sighting range 


A diagram of a spiral search is given in Figure l. 


°°. 
> 


This is essentially the same as with the additional condi- 


tion that the original target position is always investigated. 


a: 
3 
The aircraft search randomly in the circle of radius (v, en a 
max 
or the circle within which the target must be. 
o4 e 
i? 
The aircraft search uniformly over wa and (0, 27) 
The evaSive Strategies are 
Bei: 
The target chooses velocity and direction to counter oé 1 accord- 


ing to a game theory criteria. 


Po 


The target chooses velocity and direction to counter X , accord- 


ing to game theory criteria. 
CG 3: 
The target chooses speed uniformly on (0, ve ) and direction 


max 
uniformly on (0, 27). 


Py: 
The target chooses direction uniformly on (0, 27~-). The target 
proceeds radially outward so as to distribute himself uniformly over the 


geographic area bounded by r = ve colers 
max 
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Figure 1 


Diagram of Spiral Search 


Ly, 


3.3 /Spiralesearch. 

It is possible to demonstrate that a spiral search for a target 
moving radially outward from a fixed point is the most efficient way 
of searching over various target velocities, if the target direction is 
uniformly distributed on (0, 2). This may be intuitively justified as 
follows: Suppose a searcher has a given fixed range. He desires a 
search pattern which will minimize his path length while not Searching 
over any velocity range more than once. For a single circuit, i.e., a 
search of 9 from 0 to 27", the target velocities depend upon search width 
W and the initial radius of the search, To: If the velocity range is to 
remain constant throughout the circuit, the search radius muSt move out- 


ward with the velocity which would have put the target at r., which is 


0 
Lo/ty- If the searcher moves outward in such a manner, and increases the 
radial component of his velocity with each subsequent circuit such that 
r, = tT -1 ~ (v4 vA ae t, he will not re-search any previously searched 
velocity. This pattern may be flown in several ways. The outward radial 
component of velocity may be chosen to remain uniform throughout the 
search. Another method is to vary the radial component after the comple- 
tion of each circuit, thereby flying a series of spirals joined by short 
cross-legs. While this is slightly less efficient than a continuous 
Spiral, it has the advantage that it searches over a continuous target 
velocity range. This makes it more compatible with an analysis of the 
optimum velocity choices of the target, and for this reason this method 
is used. 

The equation for a single circuit is derived as follows: 


The search aircraft flies an increment of the circuit of length4 S 


in timeA t. During this time the target moves radially outward with 


l7a 


velocity Vi. 








r= vire t = Ve 2 
Sigs 
“ee 
or for a complete circuit, Ge= Si. 
Vv 
a 


To derive S, take 6 to be in radians. 





























since t= S 
Vv 
a 
hate 
In § + Cy = v, Q@ + Co 


(1) dS = rdeé 
4 Ve 
then since r = Ss 
Vv 
a 
Vv 
ds = —* sae 
Vv 
a 
or 
ee: 
95 =_t_ ag 
S Ve 
Integrating, 
Vv 
ae =f t a = 
S Vn 
or equivalently, 
Mio 
Tags 
= a + ¢C 
(2) S=cj,e 7 
: dS 
Since (1) Ee 6 
Vv 
t 
Le a 8 
(3) ae c,e 2) 
Mia 
Imposing boundary conditions, 
ee a a7 
os — 
S = 0 
vai 
C) = ig from (3) 
Ve oO 
Vv 
Co = = 2 r. .frene@oe 
we O 
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The final forms of the expressions are 





Vv 
Vv — @ 
a Va 
S = re (Ce - 1 ) 
Vet O 
Vv 
—~ @ 
r=r.e # 
O Vt 
so 
Vv 


To get these expressions in a more manageable form, express e 


asean Lntlnlee series. 








ye 
Vv 
“28 y t 9¢ 
ener = 1 + re + re + 
Vv 
a a1 
For Vv, >> Mee 
Vv 
7 ° “ 
e'a AS 1 + S) 
Vv 


Since Mies is usually over 300 and vA is usually less than 30, this 


is a valid approximation. 











Then 

Vv Vv 

3 = -Saeeeeais —— 9-1) = 6r, 
te 2 a 

and 
Vv 

r=r.(t1+ E 8 ) 

O an 


At the end of each circuit, the searcher must position himself to 
search over a new velocity increment. If the detection range either side 
of the flight path is W, the minimum velocity unsearched is (ro + W)/ ty: 
which is that velocity which would put the target at Yo + W at the start 
of the search. Since the search width is 2W, a radius of Ty + 2W will 


search velocities from (rp + W) ve ty to (To + 3W) / top, if the velocity 


is assumed to be constant across the search width. This is not strictly 


1 


To +W TO + 2w To + 3W 


true, since ———— ¢ ————— < ; however since this dis- 
to t ie 
0 0 

tance is small it will be ignored. For example, if W = 1, t, = 2, then 


0 


W/to is .5, which is the maximum deviation from the velocity at the cen- 
ter of the search strip. As the time increases, this deviation becomes 
smaller. 

If the target were at rp + 2W at to; it would move radially outward 
during the time the first circuit was being flown. The time required to 


search the first velocity increment 1s 





Then the initial radius for the following circuit is 


to ar 


) 


ry = To + 2W + 





The distance from the position at the end of the first circuit to the 


beginning point of the next circuit is 


2 mil) aw To 
Ee oy Ot + eee 0 
Yao tov, 
ie Ty = 100, W= 1, ls 300, and Ge = 22 
8 +100 
= 2 eee a Oe 
‘i * 37300 


Since these cross-legs are small compared to total aircraft range, they 
will be ignored in further computations. 
Then to compute any r, it 1S necessary to take its position at to 


and to add to it the distance traveled while the aircraft was searching 


previous circuits. n=l : 
Cn + 2nwW aa 1 
Pr, = 0 + anw + ( 0 ie) 
ty Vv 
a 
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since 


s =2W,r ; 
n n 
n-l 
ro + 2nW ae 
S =2 |r. + 2nw + ( — 
n O t Vv 
O a 
n-1l 
Because of the term Pa Si a recurSive relation exists. 
i=0 
A solution for N =n is possible if the total path flown is re- 


max 


stricted by the aircraft range, or 


N 
Re 2) Be 
i=0 


However, because of the complexity of the expression an analytic solution 
is very difficult and an iterative method of solution for the maximum 
velocity searched as a function of To is laborious. For these reasons a 
computer program was used to solve the expression. Thus a method exists 


for solving for the range of target velocities searched for any r., with 


0 


input parameters W, R, v_,; Vv, and t,. This solution is necessarily 
e max 0 

discrete, since the r9o's are discrete; however, any degree of accuracy 
desired is obtainable. 
3.4 Game Theory Evaluations of the Proposed Strategies. 
3.4.1 Discussion of Strategy + 

Using the previously derived expression for the range of velocities 
searched over for any Lo» the minimum and maximum target velocities may 


be obtained. The probability of detection for any given initial Search 


* 
radius, say ty? 1S 


x & _ _* 
PL vn (rQ) ¢ ak Vmax (rQ) / oo cf 


rag 


Then 


MCA, , Py) = P (D) 


=i Pixs Ge) < v,< ae cS) | ro rt) («, = r* ) 

The lower the value of rH the greater is the range of target vel- 

ocities searched. For the last circuit, generally not all values of 8 

are searched since the circuit must terminate when maximum range is reached. 

The target direction of travel is distributed uniformaly on (0, 2 ), which 
means the probability of detection on the last circuit is the percentage 


of the circumference Searched. This probability is denoted 8,- The den- 


Sity function for P(D/ro = aA) is 


eS = 





The graph of this density function is depicted in Figure 2. 

Denote the range of velocities searched by Vis The payoff function 
: a ae _ _* » = 
is M (vp> v> = PCD / ae and ve Me): 

The distribution functions for ve and Vv, are F(v_) and Gv, )-. These 
are not fixed distributions, but are any possible distributions. The ex- 


pected value of the game for two distributions F and G is denoted E(F,G). 


uae Me 
max max 
EC(F,G) = { MCv v,) dF(v_) dG(v,) 
0 0 
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P(D) 


oO Ose SE | 


G, 
O = (ot eNwoy C+ anwew 
oe ~ at Aare 
c } 
oO 
Figure 2 


Conditional Density Function of the Probability of Detection 
Le M(v ~~ is continuous then a solution is guaranteed. However, in 
the case of spiral search M is not continuous, and therefore no solution 
is guaranteed. However, an examination of the proof of this theorem in 
4 indicates that the requirement of a continuous payoff function may be 
too restrictive. The assumption of continuity 1s necessary in the proof 
to guarantee the existance of the Stieltjes integral of 


f mex, y) dF(x) 


and 


J ce, yeaa cy) 


Zs 


4s 


and to guarantee the existance of certain minima and maxima over the do- 
main of the payoff function. It has not been determined if the existence 
of a solution to the problem under discussion exists. 

The existence of a game theory solution has been presumed and den- 
sity functions for the searcher and target derived. Due to the complex- 
ity of the analytic expression giving the range of target velocities 
searched for a given To. an approximation to the solution using discrete 
values of Co and v, has been used. 


t 
Le ro and v, are allowed to take on discrete values and the probabil- 
ity of detection derived for each combination of these discrete values, a 
matrix approximating the continuous-domained payoff function may be ob- 
tained. The degree of accuracy is dependent on the fineness of the incre- 
ments of a and ra? for purposes of the present investigation increments 


of one unit of speed, 1.e., one knot or one mile per hour were used. The 


resultant matrix 1s m x m+l, where m = Mae . There are several meth- 
max 


ods to obtain the game theory solution of this matrix, but the precise 
methods of solution have the detriment of being extremely complex for a 
matrix of any size. For this reason an algorithm which gives an approx- 
imate solution was used. The algorithm and the computer program used to 
perform the calculations are discussed in Appendix I. An extenSive dis- 
cussion with proof of convergence is given in Karlin [2]. The relevant 
properties of the method are: 

(1) An approximate strategy for A and B may be found. 

(2) Upper and lower bounds on the value of the game are generated. 
The accuracy of the solution depends on the number of iterations per- 
formed. Two thousand iterations are considered to give a usable approx- 


imation for the matrix under discussion. As an accuracy check one program 
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was run which performed twenty thousand iterations. The upper and lower 
values of the game, denoted V and V, were .1254 and .0845 for 2000 iter- 
ations and .1153 and .1017 for 20,000 iterations. The total computer 
time time for the more accurate method is four minutes and seven seconds. 
With 2,000 iterations ten sets of input parameters may be run in three 
minutes and fifty seconds. Using this program optimal strategies for A 
and B were found, with A restricted to using only spiral search. The 


strategies for searcher and target for varying v and t. are tabulated 


Cnax 0 


in Appendix II. 
3.4.2. Discussion of Strategy of 5 

All of the theory and general discussion relating toA, is appli- 
cable toe 9° The same computer program with appropriate modifications 
was used to solve the payoff matrix. The range of the search aircraft 
was diminished by the radius of initial search to compensate for search- 


ing v, = 0 on each search pattern, which had little effect on the out- 


t 
come. This strategy forced B to always move from the position of initial 
detection, which to some extent alleviates the assumption that the time 
the target begins its movement is a deterministic rather than a probabil- 
istic quantity. The reason for this assumption is that there is a lower 
bound on the length of time the site takes to disassemble components and 
makes ready to move. If B is forced to move, the longer he is in motion 
the more difficult is the problem of detection. If B knows the initial 
position will always be searched, he is more likely to make an effort to 
become highly efficient in preparing to move, and therefore the range of 
times at which he moves cluster more closely around the minimum move time. 


The solution generated by the computer for this strategy indicates that B 


will move at a very low speed in about the same proportion that he did not 
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move when A plays strategy a 1° 
3.4.3. Discussion of Strategy OX .. 
This strategy requires A to search randomly in a circle of radius 


(tg ° Bi ). It has the advantage of not requiring any expenditure of 
max 


effort on navigation, and in fact is often used in practice. It is an 
inefficient search method and is included largely as a basis of compari- 
son for more effective strategies. The area in which it searches is not 
constant if the target is moving radially outward, and the probability 
of detection is less than is given in the following formula. However, 
since the probability of detection is low for this method, no accuracy 
was lost in the general conclusions. 


The probability of detection is given by 
-2WR 
A 
P(D) = l-e (3) 
The derivation of this formula is given by Koopman. fa] 
3.4.4. Discussion of Strategy *, . 
In this strategy A searches all velocities of B with equal likeli- 


hood, or ae is distributed uniformaly on (0, ve ). There is a problem 
max 


in this strategy arising from the fact that searching some velocities 
requires less distance traveled than others. The strategy is therefore 
formulated as the sequential choice of velocities to be searched with 
each choice being uniformly distributed on the target velocity range. 
The expected probability of detection based on the expected range of 
velocities searched is used in the payoff matrix of the various strate- 
gies. 

3.4.5. Discussion of Strategies PB , and f 5° 


These strategies are defined in the course of deriving strategies 
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HX and &,. As with &, and & 


1 D discrete values are used to approximate 


Pigs 


the continuous strategy. 
3.4.6. Discussion of Strategy p 3° 


B moves at time t, with a velocity chosen uniformly on (0, es ) 


max 
and continues to move with this velocity throughout the time period. The 


0 


direction of travel is chosen uniformly on (0, 27). This distribution 

is also obtained if the target moves with velocity wa and stops at 
max 

times uniformly distributed on (0, ty) where O denotes the minimum time 

at which the target may move. 


3.4.7. Discussion of StrategyB wn 


In this strategy player B distributes himself uniformly over the geo- 


graphical area of radius a - to- The probability of detection is 
max 
a 
PCD — ee 
(D) rn 
where 
a= Wr? 


The distribution function in terms of radius, r, is 


2 2 2, 
F(r) = P(D) ='|——— 8 —— 
t 
tmax 9 cae. 0 
The density function is 
f(r) = dF(r) = — dr 
(v ty) 
tmax 9 
Since r = tov» 
dr = tydv 


The density function in terms of v is 


Ge. 2tgvt dv 7 Oy ape 
=< Seek a Aliens) i 
V t V 
[Eo aagl t 


m 


ah 


andwene Gistribwmtion is 


if 2 
r V 
| 
P< a) = i ae 
0 ee 


Since strategies A | and sa have been approximated using discrete vel- 
ocity increments this strategy may be defined in velocity increments using 


the distribution function, which gives an expression compatible with | 


and & 5° 


28 


IV. Evaluation of Strategies. 
4.1. Method of Evaluation of Strategies 

The pairwise comparison of all combinations of strategies considered 
may be made and the results tabulated in a matrix of payoffs. Examina- 
tion of this matrix then reveals the optimum or combination of optimum 


Strategies. The elements of the matrix are 


2, = Piple =o, 9 B -B J 
The general expression for the probability of detection for strategies 


fone and & , is 


1 
eee Kore 
P(D) = S. P(rg = cr.) P(D / To = = =F P(ro = LANG 
i=1 af 


Ce Se “t Vmax (tr) 


The comparison of & l and B ? and of of D and 3 9 have been made in the 


solution of the game. The solution for elements Ar and ay, may be made 


using the general formula directly. To solve for elements 443? Ayr 259) 


aoy the following modification of the basic formula is used: 


sane Vmax ‘F4) 
P(D) = oa Gas cr.) f. £(v,) HOt (Se and Pi 
ie | : 
* Vmin (r5) 
where 
Vmax v2 - ve 
Vt 
Ymnin max 
and 
“a x 
a V - Vv: 
max min 
i BCS aor aaeall for f 3 
Vow max 
min 


Zo 


For & , the formula for the probability of detection of a uniformly dis- 
tributed target is used although this is not true for strategies ? 

F , and p 3° However, the low probability of detection for random search 
indicates that more precise computation is not necessary. 

For b mn there is a difficulty in that path length for searching dif- 
ferent velocities is not constant. This has been handled by using the 
expected value of the path length for a single circuit. The expected 
value of va is calculated for each strategy of B and this is used to get 
the expected number of circuits. This is then used to calculate the prob- 


ability of detection. 


The expected target velocity is 


Nit 
max 
HCL, Wage ORES. ) 
L L 
0 
Then 
E [path length] = E(u) = 2 Ft, E(v,) 
Thus 


E {number of circuits} = E[c] ~ oy 
While it would be possible to solve for the exact probability of detec- 
tion of & 4 for hs 1 and p 5 this was not done since rough calculations 
indicated that such a value would be very low. The expected number of 


circuits may be used to solve for an directly by 


aw 


Vee 
max 


since the target speed is distributed uniformaly on (0, v The ex- 


Ie ie 
max 


pected probability of detection for a), has been calculated to be 
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o |* 


| ey 


Vv 





t 
max 


for a Single circuit. Therefore ay), was taken to be 


W 
272 to 
es: “GRKG)) 
Vv 
t 
max 
The payoff matrix for parameters R = 500, ty = ya (i vy = 30 
max 
and vo = 300 is 
1 2 3 4 
; 108 103 . 180 Oval 
5 -103 101 174 -074 
4 O41 041 041 O41 
F 068 064 = O97 -054 


This matrix has a saddle point at a Therefore pure strategies exist 


24° 
for players A and B. It appears that searching by strategies 6 | and 

X, is highly desirable for player A, since the values of a1, and a0; 

are very close and a deviation by B from strategy B 4 greatly imeneuaae 
the probability of detection. If B distributes himself uniformly over 
his range of velocities then A increases the probability of detection 

by more than a factor of two if strategies 1 and &. are used. To 
determine whether the existence of a saddle point in this case was unique, 
values of the elements of the first two rows of the matrix were calcu- 


lated for six other sets of input parameters and in each case a saddle 


point existed at either ayy OF ay, - 
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V. Conclusions. 
5.1. Evaluations of Strategies 
The low value of strategies 3 and ng the random and uniform on 


v, searches, indicate that their use in an actual situation would not be 


t 
optimal. However, the use of a simple strategy that corrects some of the 
deficiencies of these strategies may be used. If an area search is used 


which does not re-search any area, the probability of detection will be 


approximately 


This is very slightly less than the probability of detection using strat- 


egies ©<_ and x, the spiral searches, and is very much simpler to apply. 


1 
It must be concluded that a uniform geographic distribution of the 
target is optimal in the set of strategies discussed. Since no counter 
to this strategy seems to exist, this may well be an optimal solution 
over all strategies available to player B. If the equation for probabil- 


ity of detection used for non-overlapping flight path over a circle of 


radius Ue : to is used, a very close approximation to the probability 
max 


of detection given by the saddle point is obtained. This indicates that 
strategies “ , and %. are optimal in the set of strategies studied by 
virtue of their characteristic of non-overlapping flight paths rather 
than because they derive an advantage because of some other characteris- 
Gicseelt this 1s true and if Beplays strategy B Lb? uniform geographical 
distribution, the effort needed to apply strategies o 1 and ce 7 seems 
unjustified. However, the distinction between uniform geographical dis- 
tribution and uniform velocity distribution is not obvious, and conceiv- 


ably a naive player B could choose to play strategy f 3 and distribute 
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j 
if 
' 


uniformly on (0, t, ). Sinee using strategies aa or =e in this case 


max 


would considerably increase the probability of detecting the target, the 
anticipated benefits might outweigh the increased effort. This might be 
augmented by an attempt to estimate the type of strategy applied by the 
evader by goodness-of-fit tests, to determine the probability that B is 
using some strategy other than 4, and to use a counter-strategy which 


takes maximum advantage of B's mistakes. 
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APPENDIX If 
Computer Program for Solving for Strategies 1 and 2 
This program has two basic parts. The first develops the range of 
target velocities searched for various values of To using a spiral search, 
and generates a payoff matrix using increments of one unit of speed be- 
tween successive values of rpg and v,.- The second part of the program 
solves the matrix using a method for approximating the solution of matrix 


games. 


The inputs are: 


RO initial radius searched 
VA search aircraft speed 
TO length of time the target has been in motion when 


search begins 
RNG search aircraft range 


DELR amount the initial radius is incremented for each 
successive search 


W search width 
TMAX maximum target speed 


The outputs are: 


VMIN minimum target speed searched for a given RO 

VMAX maximum target speed searched for a given RO 

THETA percent of circle searched on last circuit 

RO radius corresponding to VMIN, VMAX and THETA on same 


line of printout 
A matrix of detection probabilities 


X(i) approximation of optimum percent of time the search 
should start at RO equal to i- to + W 


Y(i) optimum percent of time the target should choose 


Vis = 1-1 
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Vil upper value of game 
v2 lower value of game 


Internal symbols are: 


SUM total distance flown 
S circumference of present circuit 
SLAST percent of circumference searched on last circuit 
VI Cj) target velocity equal to (j=1) 
V 
Cmax 
SUMA (TA) ra aij EO Toe] -> SUMB(j) is maximum over j 
i=1 
V 
tmax *} 
S UMB (IB) a. ; fog 1. -3 SUMACi) is minimum over i 
b | 
j= 
NA minimum SUMB (3) 
NB maximum SUMA (i) 


YB (IE) number of times the target has picked velocity YB 
(IE - 1) 


XA (IX) number of times the search aircraft has picked initial 
velocity searched XA= IX up to present iteration. 


This program iS written in Fortran 63 for use on the CDC 1604 com- 
puter. The first part of the program uses the equation developed in 
Section 3.3 to generate minimum and maximum target velocities searched 
for various values of RO. The initial RO is an input to the program, as 
is the increment of RO, designated DELR. For purposes of this investi- 
gation increments of one unit of speed were desired, and therefore the 
Statement DELR = TO has been inserted. RO normally begins at W, which 
assures search of the initial target position. Each successive RO de- 
fines a row of the payoff matrix. After VMIN and VMAX are found, the 


Spectrum of target velocities are subjected to IF tests to find if they 
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fall within or without the detected velocities. They are assigned prob- 
abilities of detection of O or 1. Then the element of the row correspon- 
ding to VMAX is assigned a probability of detection THETA. This process 
is continued until RO = VIMAX - TO. At this point the matrix is com- 
plete. 

Part II of the program obtains an approximate game theory solution 
for t is matrix. Obtaining an exact solution for a large matrix is diffi- 
cult unless a saddle point exists, which requires that 


max min aig = min max ai 
cl a j a 
But for each j in the matrix there exists an element Bee such that 
aij = O. For each i there exists an an such that 2 QO. Therefore 
no saddle point exists. 
The exact methods of computation involve examination of all adjoints 
of the submatrices of the matrix. The generation of the adjoints re- 


quires evaluation of the inverse and the determinant of the submatrix. 


30 
The number of submatrices of a 30 x 30 matrix is ; o (3) EPyn There- 
n=1 


fore a method approximating the actual solution has been used. This uses 
an iterative process in which player A chooses his strategy by looking 

at the column vector which is the sum of all columns chosen by B, and 
playing the row corresponding to the maximum element in this vector. B 
looks at a row vector composed of the sum of all rows chosen by A and 
chooses the minimum element for his next play. The number of times each 
strategy is chosen divided by the total number is an approximation of 
the optimal strategies. The minimum element in the row vector and the 
maximum element in the column vector divided by the number of iterations 


give lower and upper bounds on the value of the game. 
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The row vector of sums is SUMB in the program, and the column vector 
of sums is SUMA. The number of times each row and column has been chosen 
are XA and YB. Statements 70 through 80 choose the maximum over SUMA. 
Statements 74 through 77 choose the minimum over SUMB. 

The program as presented generates the solution for strategy & ,. 

To obtain the matrix and solution for strategy % j? Statement 97 and the 
preceding statement may be deleted, along with the fourth statement after 


statement 1/7. 
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